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The persistence properties of a set of random walkers obeying the A + B ^ 9 reaction, with equal 
initial density of particles and homogeneous initial conditions, is studied using two definitions of 
persistence. The probability, P{t), that an annihilation process has not occurred at a given site has 
the asymptotic form P{t) ~ const + t~^, where 9 is the persistence exponent ("type I persistence"). 
We argue that, for a density of particles p ^ 1, this non-trivial exponent is identical to that governing 
the persistence properties of the one-dimensional diffusion equation, dt(t> ~ Oxxtj), where 9 ~ 0.1207 
0. In the case of an initially low density, po <C 1, we find S ~ 1/4 asymptotically. The probability 
that a site remains unvisited by any random walker ("type II persistence") is also investigated and 
found to decay with a stretched exponential form, P{t) ~ exp(— const x Po^'^t^^'^), provided po ^ 1- 
A heuristic argument for this behavior, based on an exactly solvable toy model, is presented. 



I. INTRODUCTION 



Diffusion-limited reactions are used to model a wide 
range of phenomena in physics, chemistry and biology, 
and continue to stimulate current research. Reaction- 
diffusion processes have been applied to studies such as 
stochastic spin- flip dynamics , exciton-exciton dynam- 
ics in tetra-methylammonium manganese trichloride 
the kinetics of bipolymerisation Q , reptation of DNA in 
gels 1^, interface growth diffusion of zeolites [Q and 
other phenomena such as competing species in biology, 
self-organised criticality, pattern formation and dynamic 
phase transitions. 

Much of the effort to date has focused on reactions of 
the form A + A (d and A + B 0, with a variety of 
boundary and initial conditions |Q . It is well appreciated 
that within the context of these models there exists an 
upper critical dimension dc, below which spatial fluctu- 
ations in the initial distribution of the reactants play a 
significant role in the evolution of the density of the par- 
ticles. This dependence on the microscopic fluctuations 
invalidates traditional approaches such as the mean-field 
approximation. Attempts to understand the role played 
by fluctuations have involved numerous techniques, in- 
cluding Smoluchowski-type approximations Q and fleld- 
theoretic methods p-pT|. In this paper we set out to 
elucidate the effects of these fluctuations on persistence 
within the context of the A + B ^ (l) model. 

Persistence phenomena have received considerable at- 
tention in recent years . Theoretical and com- 
putational studies include spin systems in one and 
higher jl^ dimensions, diffusion fields |lHl^, fluctuating 
interfaces [0 and phase ordering dynamics Ex- 
perimental studies include the coarsening dynamics of 
breath figures ll9|] , soap froths and twisted nematic 
liquid crystals pi| . Persistence in nonequilibrium crit- 
ical phenomena has also been studied in the context of 
the global order parameter M (t) regarded as a stochastic 
Reaction-diffusion models offer much scope 



process |22 



for the study of persistence and have already contributed 



significantly to our understanding in this area |23|| . 

The definition of persistence is as follows. Let ip{x, t) 
be a stochastic variable fluctuating in space and time ac- 
cording to some dynamics. The persistence probability is 
simply the probability P{t) that at a fixed point in space, 
the quantity sgn [(/){x, t) — {4>(x, t))] does not change up to 
time t. In many systems of physical interest a power law 
decay, P{t) ~ , is observed, where 9 is the persistence 
exponent and is, in general, nontrivial. The nontriviality 
of 6 emerges as a consequence of the coupling of the field 
(t){x, t) to its neighbours, since such coupling implies that 
the stochastic process at a fixed point in space and time 
is non-Markovian. 

Two distinct types of persistence emerge naturally in 
the study of the A + B % reaction-diffusion model. 
Consider a nonequilibrium field (j){x, t) which takes values 
at each lattice site x. The field evolves in time t through 
interactions with its neighbours. In type I persistence, 
the field (j){x, t) changes its sign whenever an event oc- 
curs at the lattice site x at time t, where an event is 
defined to be the reaction process A + B — > 0. Type II 
persistence satisfies the conventional definition and the 
field (j){x, t) changes sign when the lattice site x is visited 
by either an A or particle at time t. The persistence 
probability at time t is defined as the fraction of sites in 
which the stochastic field 0(x, t) did not change its value 
in the time interval [0,t]. Our analysis suggests, in the 
case of type I persistence, that when the initial density 
of particles is high (po S> 1, where po is the density of 
either A ot B particles), the value of the persistence ex- 
ponent is identical to that which emerges in the study 
of the one-dimensional diffusion equation dt4> = dxx't'^ ^-s 
long as the running density, p{t), satisfies p(t) ^ 1. This 
yields 9 ~ 0.1207 When, however, the initial density 
of particles is low (po ^ 1), we find ~ 1/4, in agreement 
with a simple heuristic argument based on the decay of 
the density. In the case of type II persistence we show 
that the persistence decays with the "stretched exponen- 
tial" form P[t) ~ exp(— const x pj^^t^/**), and we provide 
a heuristic derivation of this based on the behavior of a 
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toy model. 

The paper is organised as follows. In section I we in- 
troduce the model and calculate the evolution of the par- 
ticle density for all po, including po ^ 1. In section II 
we present our results for type I persistence. In section 
III we introduce a generalised toy model of noninteract- 
ing diffusing particles and illuminate, through a special 
case, the type II persistence properties of the A + B ^ (d 
model. All our predictions are tested by extensive nu- 
merical simulations. 



II. PARTICLE DENSITY IN THE A 
MODEL 



B 



We are concerned with the following model. Consider 
the A + B ^ (d reaction involving two types of parti- 
cle, both executing diffusive random walks. The par- 
ticles move on a one-dimensional lattice with periodic 
boundary conditions, and react upon contact to form 
an inert particle. At t = 0, exactly equal numbers, 
A^a(O) = Nb{0), of A particles and B particles are ran- 
domly distributed on the lattice. This is done by ran- 
domly assigning one half of the lattice sites to A parti- 
cles and the other half to _B-particles. The A-sites and 
B-sites are then randomly filled with A and B particles, 
i.e. each A(i?)-particle occupies each of the A(i?)-sites 
with equal probability such that, at large scales, both 
densities pa{Q) and pb{Q) are initially homogeneous. We 
define po — Pa{0) — Pb{0) = A'^A(0)/i where L is the size 
of the lattice. The two species are also given the same 
diffusion constant Da = Db = = 1/2. Our model 
then evolves permitting multiple occupancy of sites, but 
we impose an instantaneous reaction so that each lattice 
site contains only one type of particle. 

The mean-field approach yields an asymptotic decay 
of the particle density according to pA{t) = PB{t) 
const X with an amplitude independent of the ini- 
tial density. However, in low enough spatial dimension, 
d < dc = 4, the dominating process asymptotically is the 
diffusive decay of the fluctuations in the initial condi- 
tions. This leads to anomalous kinetics and it was shown 
by Toussaint and Wilczek that 



{pAit)) = {pB{t)) 



7,1/2 (8^)<i/4 



(Dt) 



-d/4 



(1) 



This result has been confirmed using field-theoretic meth- 
ods for 2 < d < 4 by Lee and Cardy The exponent 
—d/A has also been rigorously confirmed by Bramson and 
Lebowitz pst . The numerical simulations which have 
been performed in one |p^ , p6|| , two and three di- 

mensions [ p8[ are also in good agreement with the an- 
alytical predictions of the decay exponent. However, it 
was noted by Lee and Cardy that, although in one 
dimension reasonable agreement with the dependence on 
the initial density has been found, in higher dimensions 
the \/p{0) amplitude dependence has not been observed. 



They suggest that, in the one-dimensional simulations, 
the initial average occupation number per site was kept 
low, whereas for the higher dimensional simulations it 
was necessary to start with a nearly full lattice in order 
to reach the asymptotic regime and therefore that Eq. 
(|l|) might not be a universal result, but rather a limit 
for small initial density. We believe this to be the case. 
We focus on the particle density in one dimension. In 
our study of type I persistence we are particularly inter- 
ested in the limit of high initial particle density, po 3> 1, 
where the problem should map onto simple diffusion. As 
a precursor to looking at the persistence, therefore, we 
first consider how the particle density decays with time 
in this limit. 

In the high-density limit, where there are many par- 
ticles per site, one can neglect dynamical fiuctuations in 
the density. If Nn{t) is the number of particles at site n 
at time t, the dynamics of Nn(t) is governed by diffusion: 



Nn{t) = D(Nn+i{t) - 2N.„{t) + Nn-i{t)) 



(2) 



Here we have adopted the convention that Nn > means 
that site n is occupied by iV„ A-particles, while iV„ < 
means that it is occupied by —Nn B-particles. Then the 
annihilation process is automatically built into the diffu- 
sion equation (|^). Introducing a discrete Fourier trans- 
form and taking the limit L oo, gives, 

iV„(t) = Va^™(0) / ^exp[-2I?(l-cosfc)i 

+ik(n — m)]. (3) 

For t 3> 1 we are justified in using the approximation 
I — cos/c ~ fc^/2 in the integrand, and extending the 
limits on the fc-integral to ±oo. This gives 



in — to)^ 



(4) 



For i — > cx), the Gaussian kernel becomes broad, and 
Nn(t) becomes, asymptotically, a Gaussian random vari- 
able (where the randomness comes from the initial con- 
ditions). It's mean is clearly zero, since (iVm(O)) = by 
symmetry, so its distribution is completely specified by 
its variance. This is independent of n by translational 
invariance, so we drop the subscript and write (again for 
t oo). 



E 

ni.m' 



X exp 



(A^„(0)A^„KO)) 
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ADt 



(5) 



Using {N,n{0)N,n'0) = {N^^{0)) Smm' wc obtain, asymp- 
totically. 



^/SirDt ' 



(6) 
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At time t = the number of particles of a given type 
on a given site has a Poisson distribution with mean 2po, 
where recall that po is half the total density. Hence the 
variance is given by {N^{0)) = {\N{0)\'^) = (2po)^ + 2po, 
giving 



{N\t)) 



a\t) 



(7) 



where a is the standard deviation, and the prob- 
ability distribution of N{t) is simply P[N{t)] = 
[l/V^h^) exp(-7VV2a2). 

The total particle density, recalling that, with our con- 
vention, i?-particles have been assigned negative values 
of Nn{t), is given by 



dNN exp 



2^ 



a. (8) 



Substituting for a from Eq.(0) yields, finally, 



{\Nm = 



(9) 



The individual particle densities, therefore, decay asymp- 
totically according to, 

= ip.it)) ^ M±0^(i^i)-V4. (10) 

It is interesting that in the low density limit pp ^ 1, 
we recover the result of Toussaint and Wilczek [E4| , 



{pA{t)) = {pB{t)) 



1/2 

Po 
(27r)3/4 



(11) 



However, in the high density limit po ^ 1 we obtain 
the same decay exponent —1/4 as for low density but a 
different amplitude. The amplitude now scales as po in 



contrast to the low density limit where it scales like Pq 



1/2. 



{pA{t)) = {pB{t)) 



PO 



(27r3)i/ 



-1/4 



Po > 1- 
(12) 



Given the correctness of Eq.(pX|) in both low and high 
density limits, we expect it to be a good approximation 
across the whole range of densities. 

For convenience, let us write Eq.(p^) as (p^(t)) — 
(psit)) — K pQ{Dt)~'^ . The analytical values of K and 
7 are then K = 0.35635 . . . and 7 — 1/4. Numerical 
evidence has already substantiated Eq.(|ll]) |^,^. Be- 
low, we present numerical support for our high density 
calculation. 



2 - 



p„=500 

Po=250 

Po=50 







10 



15 



Int 



FIG. 1. Log-log plot of the particle density as a function of 
time for po = 50, 250 and 500. 



Po 




K 


50 
250 
500 


0.2576 ±0.0005 
0.2524 ±0.0002 
0.2548 ±0.0005 


0.3599 ±0.0007 
0.3568 ±0.0003 
0.3669 ±0.0008 



TABLE 1. Numerical values of the decay exponent 7 and 
the constant K where {pA{t)) = {pB(t)) ~ Kpo{Dt) ' for 
Po ~ 50, 250 and 500. The numerical values of K were ob- 
tained by plotting {pA{t)) vs. t ' and evaluating Knum from 
the gradient of the curve. The analytical values are 7 = 1/4 
and K = 0.35645. 

In conducting numerical simulations of this nature one 
must naturally consider how to increment time. A phys- 
ically realistic and computationally simple method is to 
increment dt by (1/current number of particles). How- 
ever, the expediency of this technique clearly does not 
extend to systems with high densities of particles. A 
more efficient method is to allow all particles to move 
simultaneously. In our simulations, therefore, one time 
step constitutes a jump by all particles in the system to a 
nearest neighbour site with equal probability. In the low 
density limit, we have made comparisons of the two meth- 
ods of updating and find, asymptotically that the results 
are identical. Within the context of a low initial density, 
systems which permit multiple occupancy of lattice sites 
generate indistinguishable results from those which per- 
mit a maximum of one particle per site. For the purposes 
of this paper, therefore, all of the simulations are based 
on a parallel updating of the particles and permit multi- 
ple occupancy of lattice sites. Our numerical simulations 
for an initial high density of particles are performed on 
a one-dimensional lattice of size 10*. Each run is per- 
formed for 5 X 10* time steps and we average our results 
over 100 runs. The results are shown for three different 
initial densities in Figure 1. After initial transients, the 
log-log plots apparently approach straight lines in each 
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case. The values extracted for the decay exponent 7 and 
the ampUtude K are given in Table 1, where the quoted 
errors are purely statistical. We attribute the small dif- 
ferences between the measured and theoretical values to 
a failure to reach the true asymptotic regime. 



III. TYPE I PERSISTENCE 



Consider the A + B ^ % reaction-diffusion process 
on a one-dimensional continuum. The rate equations 
for the concentrations are OiNa = dxx^A — R and 
dtNB — OxxNb — R where R is the reaction rate per unit 
volume. The concentration difference, AiV = Na — Nb, 
obeys the simple diffusion equation dtcj) = d^x'P- It is 
well known that, below the critical dimension dc — 4, the 
A + B ^ % model evolves to a coarsened state in which 
the two species segregate into domains of either A oi B 
particles. The domain walls are defined by AiV = and 
their motion is clearly determined by the annihilation 
process A + B 0. The persistent sites in this model 
are therefore defined as those sites which have not been 
touched by the zeros, AiV = 0, of the diffusion field, i.e. 
those sites which have never seen an annihilation pro- 
cess. In order to overcome the discrete nature of our 
problem and model the diffusive process of the domain 
walls as described above, we allow our model to approach 
the continuum limit by starting with a very high initial 
density of particles, po 3> 1. One might think that, in 
this limit, the entire system will asymptotically become 
nonpersistent. However, we find this not to be true. The 
persistence decays according to P{t) ^ A + t^^ where A 
is a constant i.e. there are some sites which are always 
isolated from the diffusing boundaries between domains. 
We expect the offset A to vanish as po 00: clearly, in 
the real, continuous diffusion problem P{t) ^ t^^, with 
no offset. Our hypothesis is that, in the limit of a large 
initial density of particles po ^ 1, the offset tends to zero 
and the persistence exponent is identical to that which 
emerges from a study of the diffusion equation in one di- 
mension, 6 ~ 0.1207 [0. This should hold if the density 
during the entire run remains large enough {p{t) ^ 1) 
that the continuum description in terms of the diffusion 
equation remains valid. If, on the other hand, the sim- 
ulation enters the low density regime at late times, we 
expect the approach of the persistence to its asymptotic 
value to be described by the low-density theory. 

In the limit of low initial density, po 1, sl realistic 
representation of the diffusing boundaries according to 
the diffusion equation is not feasible, due to the discrete 
nature of the problem. In this case it is clear that the 
entire system cannot become persistent since, by defini- 
tion, the total number of reactions A + B ^ (!) L 
where L is the size of the lattice. Hence the relationship 
P{t) A + t-^ is expected, with A ^ 1 as po ^ 0. The 
large mean separation between particles implies that any 
site which experiences an annihilation process is unlikely 



to witness such an event again i.e. every reaction pro- 
cess A + B ^ ^ leads to a single, isolated nonpersistent 
site. Given that the total number of reaction processes 
is governed by the initial density of walkers, in the limit 
of sufficiently low density one can set a lower bound on 
the persistence {P{oo)) ~ A > 1 — po- The persistence 
properties are therefore determined by the decay of the 
walker density. The walkers decay according to ~ t~^^^ 
and we therefore expect that, for po 1, 9 — 1/4. 

Our numerical results are presented below. The alge- 
braic form for the persistence prevents a simple evalua- 
tion of the exponent. A three parameter nonlinear curve 
of the form A + Bt~^"' was fitted to the data to ascertain 
the value 9ni where the subscript denotes the implemen- 
tation of a nonlinear curve fitting technique. The fit was 
carried out using t as the independent variable, though 
the data are plotted against Int for clarity of presenta- 
tion. To demonstrate the validity of determining 9 in this 
manner, we also plot our data in the form {P(t)) against 
t'^"' , where the correct value of 9, in such a plot, mani- 
fests itself as a straight line graph. An unbiased method 
of determining the exponent is also presented, by numer- 
ically differentiating the data and displaying the results 
on a log-log plot. The value of 9 determined in this man- 
ner is denoted 9diff. We consider the high density limit 
first and take as our initial values po — 500 and 250. The 
simulations are performed on a one-dimensional lattice of 
size 10"* and run for 5 x lO"' time steps. We average our 
results over 100 runs. In the nonlinear fits, the regime 
i < 3 X 10^ is fitted, corresponding to Int < 10.3. 
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FIG. 2. A nonlinear fit of the form A + Bt~ is made to 
the fraction of persistent sites at time t for an initial density 
po = 500. The values of the fit parameters are A — 0.05±0.02, 
B = 0.57 ± 0.005, and e„i = 0.13 ± 0.01. 

The data are presented in Figures 2-7. Figure 2 shows 
the nonlinear fit for pQ — 500, from which the expo- 
nent 9ni = 0.13 ± 0.01 is extracted, consistent with the 



4 



diffusion value, 9 ~ 0.1207. The plot, in Figure 3, of 
(P(t)) against ^-0129 (where 0.129 is the best-fit value 
of Oni from Figure 2) yields a good straight line at large 
t. The unbiased determination of 9 from a numerical 
differentiation of the data is shown in Figure 4, where 
the logarithmic derivative of {P{t)) is plotted against 
f on a log-log plot. The resulting data should be a 
straight line of slope —9. The best-fit value of the slope 
is —9diff = — 0.12± 0.01. The slope was extracted from 
the region between the arrows, where initial transients 
have decayed but the data are not yet too noisy. The 
expected slope, —0.1207, is shown as a guide to the eye. 
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FIG. 3. The fraction of persistent sites is plotted against 
f-eni for po = 500. The value e„i = 0.129 is taken from the 
nonlinear curve fit in Fig.2. 
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FIG. 4. The fraction of persistent sites at time t is numeri- 
cally differentiated (A (Int) =0.1) and presented on a log-log 
plot for po — 500. The gradient, taken between the arrows, 
gives Odiff = 0.12 ± 0.01. The straight line has slope -0.1207. 

The equivalent data for po = 250 arc shown in Figures 
5-7. In this case one obtains 9ni = 9diff = 0.13 ± 0.01, 
again consistent with the diffusion equation result. 



0.7 



0.6 



0.5 



D. 

V 



0.4 



0.3 



0.2 



P„=250 



data points 
fitted curve 







FIG. 5. Same as Fig.2, but for po = 250. The values of 
the fit parameters are A = 0.07 ± 0.02, B = 0.55 ± 0.005, 
0„i =0.13 ±0.01. 
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FIG. 6. Same as Fig.3, but for po = 250, with Oni = 0.132. 

Notice that the values of A obtained from the non- 
linear fits are quite small - 0.05(2) and 0.07(2) for 
Po = 500 and pg = 250 respectively. An equivalent fit 
for Po = 50 gives the larger value A = 0.17(2), so the 
data are consistent with the hypothesis that ^4 — > for 
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Po — > CO, although the data are not good enough to ex- 
tract the functional dependence of A on po in this limit. 
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FIG. 7. Same as Fig.4, but for po 
ii// =0.13 ±0.01. 
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= 250. 
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the logarithmic scale for the abscissa greatly expands the 
early-time regime. The nonlinear fit was restricted to the 
range 10^ < t < 10^ (i.e. the last 90% of each run), cor- 
responding to 9.2 < Int < 11.5. Figures 9 and 12 show 
the same data plotted as {P{t)) against t~^"', and re- 
veal the expected linear behaviour at late times. Finally, 
Figures 10 and 13 give the log-log plots for the differen- 
tiated data, from which the exponent estimates 6di / / are 
obtained. The slope is measured between the arrows in 
Figs. 10 and 13, corresponding to a similar region to that 
used for the non-linear fits. 
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FIG. 9. Same as Fig.3, but for po = 0.2. e„i = 0.247. 
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FIG. 8. Same as Fig.2, but for po = 0.2. The fit parameters 
are A = 0.90 ± 0.02, B = 0.06 ± 0.005, 9ni = 0.25 ± 0.01. 



Our low density results are presented below. We 
choose as our initial densities po ~ 0.2 and 0.1. In this 
regime, our simulations are performed on a lattice of 
size 10^ and run for 10^ time steps. The system takes 
longer to enter the asymptotic regime due to the low 
initial density and this is reflected in our extended num- 
ber of time steps. We average our results over 50 runs. 
Figures 8 and 11 show the same type of non-linear fit, 
(P(t)) = ^ + S/i^"' that was used in Figures 2 and 5. 
The fit works well except at early times, and the fitted 
values are close to ^ = 1/4 as anticipated. Note that 
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FIG. 10. Same as Fig.4, but for po = 0.2. Here the straight 
line has a gradient of -1/4. The gradient of the data, measured 
between the arrows, is Odiff = 0.25 ± 0.01. 

The estimates 9ni and 9diff of the persistence expo- 
nent, for both low and high density regimes, arc sum- 
marised in Table 2. A natural consequence of the slow 
algebraic decay of the particle density is that the asymp- 



6 



totic regime, in the case po ^ 1, occurs only after many 
time steps. We expect the result ~ 1/4 for the low- 
density regime to manifest itself more clearly in systems 
run for a greater number of time steps. 
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FIG. 13. Same as Fig.lO for po = 0.1. e^ff = 0.27 ± 0.01. 
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FIG. 11. Same as Fig.2, but for po = 0.1. The fit parame- 
ters are A = 0.95 ± 0.02, B = 0.05 ± 0.005, Oni = 0.26 ± 0.01. 
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FIG. 12. Same as Fig.3, but for po = 0.1. Oni = 0.264. 
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eni 


^diff 


500 
250 
0.2 
0.1 


0.05 ±0.02 
0.07 ±0.02 
0.90 ±0.02 
0.95 ±0.02 


0.57 ±0.005 
0.55 ±0.005 
0.06 ±0.005 
0.05 ±0.005 


0.13 ±0.01 
0.13 ±0.01 
0.25 ±0.01 
0.26 ±0.01 


0.12 ±0.01 
0.13 ±0.01 
0.25 ±0.01 
0.27 ±0.01 



TABLE 2. Values of the exponent 6 in the relationship 
P{t) = A + Bt~^, where has been evaluated according to a 
nonlinear curve fitting technique and through a method 
of numerical differentiation {Odiff)- We expect, asymptoti- 
cally, e ~ 0.1207 for po > 1 and 61 = 1/4 for po < 1. 



IV. TYPE II PERSISTENCE 

In this section wc address the following question: What 
is the asymptotic probability P{t) that a given site has 
never been visited by an A or B particle, in a one di- 
mensional system of A and B random walkers subject to 
the A + B ^ ^ reaction-diffusion process? Necessarily, 
we pose this question within the context of a low initial 
density. 

Our approach involves the introduction of a toy model, 
whose predictions we will compare to simtilation results. 
The model consists of a system of nonintcracting dif- 
fusing particles, with initial density po = N{0)/L, on a 
onc-dimcnsional lattice with diffusion constant D ~ 1/2. 
Note that the definition of po in type II persistence refers 
to the total density of initial particles, in contrast to type 
I persistence, where po referred to half the total density. 

Particle decay is invoked by allowing each particle 
to vanish at every time step with some (usually time- 
dependent) probability. We define f{t) to be the fraction 
of particles remaining at time t i.e. f{t) = p(t)/pQ. We 
are then free to mimic the behaviour of any reaction- 
diffusion process by a suitable choice of f{t). 

The analysis of persistence in this toy model is straight- 
forward. We first define the probability, Q{x,t), that the 
origin has not been crossed by a given diffusing parti- 
cle with initial position x > 0. We define our lattice to 
be the interval {—L/2, L/2). Then an elementary calcu- 
lation gives Q{x,t) = eri{x/2\fDi), where we have as- 
sumed that X <^ L and » Dt, i.e. we have effectively 
taken the limit L oo in the calculation of Q{x,t). A 
simple approach in the calculation of the number of per- 
sistent sites is to consider the first-passage-time distri- 
bution Pi{t), where Pi{t)dt is the probability that the 
first crossing of the origin occurs in {t,t + dt). Then 
Pi{t)dt = Q{x,t) - Q{x,t + dt) gives Pi{t) = -dQ/dt. 
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Since the diffusing particle survives to time t with proba- 
bility f(t), the probability for the origin to be persistent 
at time t is given by, 



1 



= 1 



dt'P,{t')f{t') 




ds 



f{s) exp 



ADs 



(13) 



In the presence of many random walkers, whose ran- 
dom initial positions have a uniform distribution over 
space, the mean persistence is given by {Pit)) = 
{Q{\x\,t))^°'" where pqL = N the total number of ini- 
tial walkers and the average is over the initial position of 
a walker, assumed uniform in (— L/2,L/2). This gives 



{Pit)) = 



X exp 



For L ^ oo we obtain, 




(P(<)) = exp -2po 




5I/2 



(14) 



(15) 



A system of strictly diffusing particles with no de- 
cay is modelled by setting /(s) = 1. In this case, 
(PU)) = exp(-Ati/2) where A = Apo^JD/tt. In Fig. 
n4 we present a numerical verification of our calculation. 
The gradient of the graphs gives —A. The values of A 
extracted from the data are compared with the numerical 
predictions in Table 3. The agreement is excellent. 
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FIG. 14. Log-linear plot of the persistence for a system of 
noninteracting diffusing particles with no decay. 



Po 


Ath 


A 

^num 


0.002 
0.005 


0.00319 
0.00797 


0.00320(5) 
0.00790(5) 



TABLE 3. The persistence properties of a system of nonin- 
teracting particles with zero probability of decay is described 
by {P{t)) ~ exp(— Af^/'^). Above we provide a comparison 
of the values of A as deduced from the theory Ath and from 
numerical simulations Anum- 

Our model can also be applied to the g-state Potts 
model, which has q distinct but equivalent ordered 
phases. The A + A ^ (d model corresponds to q ~ 2, 
while A + A ^ A corresponds to g = 00. The density of 
walkers in the ID Potts model is known to decay asymp- 
totically according to p9[ , 



Pit) 



1 1 



q \/2nDt 



(16) 



We can derive a generalised value of 9{q) for our toy 
model by substituting /(s) — p{s)/po into Eg. (p^) where 
p{s) is given by Eq.(n6|). This gives. 



where 



{Pit)) ~ 



9t(9) = — 



(17) 



(18) 



and the subscript T denotes 'toy'. Despite the simplistic 
nature of our model, a power law decay for the persis- 
tence is predicted, in agreement with known results. Our 
expression for 6{q), (Eq.(^8|)) is clearly a poor approxi- 
mation to the exact expression for the ID Potts model 
obtained by Derrida et al [n4| . 



<q) 



1 



V2 



(19) 



where the subscript P denotes 'Potts'. The values of 
the persistence exponent returned by our model for the 
q = 2 and q = 00 cases are, Oriq = 2) = 0.225 and 
Oriq = 00) 0.450. The values given by Eq.(|l9|) for the 
Potts model are 9p{q = 2) — 3/8 and 9p{q = 00) — 1. 
The encouraging feature of our toy model, however, is 
that the essential dynamics of the system leading to the 
power-law decay of the persistence in Eq. ([l7| ) is correctly 
identified as the t~^^^ decay in the number of surviving 
random walkers at time t. If the number of surviving 
walkers decays as for large t, the toy model predicts, 
via Eq. dl^ ) , that the persistence will decay as a stretched 
exponential, {P{t)) - {-At^/^-°'), for a < 1/2, while 
{Pit)) will approach a non-zero constant for a > 1/2. 
The borderline case, a = 1/2, yields power-law decay of 
the persistence, as we have seen. 

We naturally now turn our attention to the A -|- B — > 
reaction process to see if the correct time dependence of 
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the persistence also emerges from our model. The den- 



sity of A and B particles in this case in given by Eq. (11) 
Therefore {p{t)) is given by, 



(Pit)) 



1/2 

Po 
(27r3)i/4 



(Dt) 



^1/4 



(20) 



where p{t) — pA{t) + Psit) and po — N{0)/L. Substi- 
tuting f(s) — p{s)/po where p{s) is given by Eq.(|20|) 
into Eq.(|l|) yields {P{t)) ~ exp{-Bt^/^) with B = 

(2'^^D/TT^y/'^pl^^. For the A + B ^ $ reaction we there- 
fore predict a stretched exponential decay of the persis- 
tence, with exponent 1/4. 
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FIG. 15. Log-linear plot of the persistence for the 
A + B ^ reaction-diffusion process. {P{t)) ~ exp{—Bt'') 
with 7 = 1/4. 

The results of our numerical simulations are pre- 
sented below. The simulations are performed on a one- 
dimensional lattice of size 10^ for 10^ time steps. We 
choose as our initial densities po — 0.02 and po — 0.04. 
A direct test of the stretched exponential prediction is 
obtained by plotting \n{P{t)) against i^/'', as in Figure 
15. The fact that the data do not clearly show the ex- 
pected straight line behaviour may indicate that the den- 
sity is not yet in the regime where it is well described by 
the 

given by Eq.([20|). Indeed, a direct study 
of the density indicates that the t^^/^ behaviour is only 
evident at the latest times reached in the simulations. 
An alternative analysis involves differentiating the data 
with respect to In t and presenting the result as a double 
logarithmic plot, as in Figures 16 and 17. Here there is 
evidence that the data approach the expected slope of 
1/4 at late times, though the noisy character of the dif- 
ferentiated data at the very latest times tends to obscure 
the asymptotic behaviour. 



T3 ^ 



CL 

V 



P„=0.02 








10 



Int 



FIG. 16. The log of the persistent fraction at time t is nu- 
merically differentiated (A(Int) — 0.1) and presented on a 
log-log plot for po ~ 0.02 The gradient of the curve gives 7. 
The straight line has a slope of 1/4. 
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A 

a: -2 
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Int 

FIG. 17. Same as Fig.l6 but for po = 0.04. 

In short, the data as presented in Figures 15,16 and 17 
are not clean enough to provide a convincing test of our 
prediction {P{t)) ~ exp(-Br) with 7 = 1/4. The slow 
approach of the particle density to its asympotic l/t-^^'^ 
behaviour is, in our view, the reason for this. A better 
test of our model, therefore, is to consider the persistence 
as a function of the actual running density of particles, 
rather than using the asymptotic behaviour of the density 
in Eq. (|l^) . Taking the log of Eq. ( p5| ) and differentiating 
with respect to Int yields. 



dint V TT 



(21) 
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since pof{t) = p{t), the running density. The left-hand 



side of Eq.(|2T]) and p{t)t^/'^ (without the factor 2yjD/Ti) 
are plotted in Figures 18 and 19 for the two densities 
studied. The two curves agree rather well over the whole 
range of Int, except at the latest times, where the date 
are very noisy. The factor in Eq.(|TJ), omitted 

from the plots, has the numerical value •y/2/7r ~ 0.8 (re- 
call D = 1/2), whereas Figures 18 and 19 suggest this 
number should be closer to unity. Given the crudeness of 
the toy model, however, the agreement between the data 
and the model is surprisingly good. 
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FIG. 18. A direct test of Eq.(|l|) for po = 0.02. 



V. SUMMARY 



In this paper we have investigated the persistence prop- 
erties of the A + B ^ % reaction-diffusion process using 
two distinct definitions of persistence. In type I persis- 
tence, we studied the fraction of sites which had never 
witnessed an annihilation process. In the high density 
limit po ^ 1 we argued that the persistence exponent 
is that which emerges from a study of the one dimen- 
sional diffusion equation, dt4i = ^xx4'^ where 9 ~ 0.1207 
[|l|, and presented data consistent with this result. In 
the low density limit po 1 we argued that the per- 
sistence properties are governed by the same exponent 
that describes the decay of the particle density, giving 
9 — 1/4. Again, the data support this prediction. In 
type II persistence, we considered the probability that a 
site remained unvisited by either an A or i? particle. Our 
approach, in this case, was to develop a toy model, which 
can be applied to any reaction-diffusion process and ex- 
presses the persistence properties in terms of the density 
of the reactants. For the A + B ^ % process the model 
predicts {P{t)) ~ exp(-Bti/4). The data are consistent 
with this result when allowance is made for the actual 
time-dependence of the particle density (rather than just 
its asymptotic form). An obvious goal for future study 
would be to to place this stretched-exponential decay on 
a firmer theoretical foundation. 
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